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Abstract. In this paper we introduce the idea of probability in the definition 
of Sequential Dynamical Systems, thus obtaining a new concept, Probabilistic 
Sequential System. The introduction of a probabilistic structure on Sequential 
Dynamical Systems is an important and interesting problem. The notion of ho- 
momorphism of our new model, is a natural extension of homomorphism of sequen- 
tial dynamical systems introduced and developed by Laubenbacher and Paregeis 
in several papers. Our model, give the possibility to describe the dynamic of the 
systems using Markov chains and all the advantage of stochastic theory. The no- 
tion of simulation is introduced using the concept of homomorphisms, as usual. 
Several examples of homomorphisms, subsystems and simulations are given. 



Introduction 

Genetic Regulatory Networks had been modeled using discrete and continuous 
mathematical models. An important contribution to the simulation science is the 
theory of sequential dynamical systems (SDS) [H, El El [TDJ [TT]. In these paper, 
the authors developed a new theory about the sequential aspect of the entities in 
a dynamical systems. In particular Laubenbacher and Pareigis created an elegant 
mathematical background of the SDS, and with it solve several aspects of the theory 
and applications. 

Probabilistic Boolean Networks (PBN) had been recently introduced [T31 HU US] , 
to model regulatory gene networks. The PBN are a generalization of the widely used 
Boolean Network model (BN) proposed by Kauffman (1969) [HE]. While the PBN 
eliminate one of the main limitations of the BN model, namely its inherent deter- 
minism, they do not provide the framework for considering the sequential behavior 
of the genes in the network, behavior observed by the biologists. 

Here, we introduce the probabilistic structure in SDS, using for each vertices of the 
support graph a set of local functions, and more than one schedule in the sequence 
of the local function selected to form the update functions, obtaining a new concept: 
probabilistic sequential dynamical system (PSS). The notion of simulation of a PSS 
is introduced in Section H] using the concept of homomorphism of PSS; and we prove 
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that the category of SDS is a full subcategory of the of the category PSS. Several 
examples of homomorphisms, subsystems and simulations are given. 

On the other hand, Deterministic Sequential Dynamical Systems have been stud- 
ied for the last few years. The introduction of a probabilistic structure on Sequential 
Dynamical Systems is an important and interesting problem. Our approach take 
into account, a number of issues that have already been recognized and solved for 
Sequential Dynamical Systems. 

1. Sequential Dynamical Systems and SDS-homomorphism 

This section is an introduction with the definitions and results of Sequential 
Dynamical System introduced by Laubenbacher and Pareigis. Here we use SDS 
over a finite field. In this paper, we denote the finite field GF(p r ) by /C, where p is 
a prime number. 

1.1. Sequential Dynamical Systems over finite fields. A Sequential Dynami- 
cal System ( SDS ) over a finite field T = (T, (/j)™ =1 , a) consists of 

(1) A finite graph F = (Vp, Ep) with Vp — {1, . . . ,n} vertices, and a set of edges 
E r C Vp x Vp. 

(2) A family of local functions j\ : JC n —>■ JC n , that is 

fi\Xlj ■ ■ ■ j X n ) v^l> • • • > "Ei—li f \ •^i+l) • • • ? *^n) 

where f(xi, . . . ,x n ) depends only of those variables which are connected to 
i in T. 

(3) A permutation a = ( «(1) ... a(n) ) in the set of vertices Vr, called an 
update schedule ( i.e. a bijective map a : Vp — > Vp). 

The global update function of the SDS is / = f a (i) o . . . o f a r n \ . The function / 
defines the dynamical behavior of the SDS and determines a finite directed graph 
with vertex set fC n and directed edges (x, f(x)), for all x G /C n , called the State 
Space of T . 

1.2. Homomorphisms of SDS . The definition of homomorphism between two 
SDS uses the fact that the vertices Vp = {1, . . . , n} of an SDS and the states /C n 
together with their evaluation map 

K n x Vp 3 (x, a) \— >< x, a >:= x a G JC a 

form a contravariant setup, so that morphisms between such structures should be 
defined contravariantly, i.e. by a pair of certain maps (<f) : T — > A, : K m — > K n ) 
or by a pair (<f) : A — > T, : K n — > K m ), with the graph A having m vertices. Here 
we use a notation slightly different the one using in [H]. 

Let F = (r, (fi : /C n -> /C n ),a) and G = (A, : /C m -> /C m ),/?) be two SDS . 
Let : A — > T be a digraph morphism, and 

(0 6 : /C^( 6 ) /C 6 , V6 G A), 
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be a family of maps in the category of Set. The map h$ is an adjoint map, because 
is defined as follows: consider the pairing 

/C n x Vp 3 (x, a) >— >< x, a >:= x a G /C a ; 

and similarly 

K m x V A 3 (x, b) i->< x, b >:= x b G fC b . 
The induced adjoint map is < >:= (pb(< x,(p(b) >) = </>&(#0(&)). Then 0, 

and (0b) induce the adjoint map : /C n — > K m defined as follows: 

(I) h 4> (x u ...,x n ) = ($ bl ),..., ?6 m (^ (6m) )). 

Then /i^ is an homomorphism of SDS if for a set of orders Tp associated to (3 in 
the connected components of A, the map holds the following condition and the 
commutative diagram. 



)C n > IC n 

| \h ( f s 

I I 

jC m ► /c m 



where n^ e ^-i( ai ) IIJ=i(0/3j) = fl'A ° #A+i ' ' ' ° 9p s and {A, A+i, • • • , A} = V*)- 
If = 0, then Id/cm o = o f a . } and the commutative diagram is the 

following. 



(Dl.b) 



K n ► /C n 

I ^ | ^ 

i i 
jC m ► JC m 

For examples and properties seepj]. It is clear that the above diagrams implies 
the following one 

f=fa 1 °— °fc n 

)C n > K n 



(D2) 



| h<j, \h ( f ) 

I I 

K m > K m 

g=g 01 o-og 0rn 

2. Probabilistic Sequential Dynamical Systems 

The following definition give us the possibility to have several update functions 
acting in a sequential manner with assigned probabilities. All these, permit us to 
study the dynamic of these systems using Markov chains and other probability tools. 
We will use the acronym PSS ( or SDS ) for plural as well as singular instances. 



Definition 2.1. A Probabilistic Sequential dynamical System (PSS) 



V s = (T,{F t }^-\(^)?=i,C) 
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over K, consists of: 

(1) a finite graph V = (Vr, Ej>) with n vertices; 

(2) a set of local functions Fi = {fn~ : /C n — > /C n |l < k < £(i)} for each vertex % 
ofT. (i. e. a bisection map ~: Vr — > {Fi\l <i <n}) (for definition of local 
function see $1.11 2)). 

(3) a family of m permutations ctj = ( <x,(l) • • • <x,(n) ) in the set of vertices 

Vy. 

(4) and a set C = {ci, . . . , c s }, of selection probabilities. 

We select one function in each set Fi, that is one for each vertices i of T, and a 
permutation a, with the order in which the vertices i are selected (similarly SDS), so 
there are m possibly different update functions / = f a {i)k aW • • • ° fa(n)k a{n) > where 
m < ml x £(1) x ... x £(n). The probabilities are assigned to the update functions, 
so there exists a subset S = {/i, . . . , f s } of update functions such that Ck = p(fk), 
1 < k < s. 

The State Space of T>s is a digraph whose vertices are the elements of KJ 1 and 
there are an arrow going from x\ to X2 if there exists a functions /, such that 
%2 = f(%i)- For each one of the selected functions in S we have an SDS inside 
the PSS, so the state space of the function / is a subdigraph of the state space of 
the PSS, so, the State Space of V s is a superposition of all inside SDS . When we 
take the whole set of update functions generated by the data, we will say that we 
have the full PSS. We denote by S the complement of set S in the set of all update 
functions, and we will call the PSS T>^ building by the same data but taking S as a 
set of update functions, one complement of the PSS Vs- All the complements have 
the same set of function but they can use different set of probabilities. 

The probability of the arrow going from x\ to X2 for example, is the sum of 
the probabilities of all functions /, such that 22 = fi^i)- The PSS represents a 
generalization of the SDS : a SDS is a PSS for which every set of local functions has 
one element, and there is only one permutation in the family of permutations. The 
update functions of the PSS have probabilities and the state space of the PSS (or 
high level digraph) is described by a transition matrix, and the dynamic is described 
by a Markov Chain. 

2.1. Example. Let 

V = (V; Fi, F 2 , F 3 ; a 1 , a 2 ; F, (C(/ f ))? =1 ), 

be the following PSS: 

1 — 3 

(1) The graph: T \ \ . 

2 

(2) If x = (xi, x 2 , x 3 ) G {0, l} 3 , then the sets of local functions are the following: 

^1 = {/n( x ) = (!, x 2 , x 3 ), /i 2 (x) = (x x + 1, x 2 , x 3 ))} 
F 2 = {/2i(x) = (x 1 ,x 1 x 2 ,x 3 )} 

F 3 = {/si(x) = (Xx,X 2 , XiX 2 ), / 32 (x) = (Xi, X 2 , XiX 2 + x 3 )} 
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(3) The schedules or permutations are ai = ( 3 2 1 ) ; a 2 = ( 1 2 3). We 
obtain the following table of functions, and we select all of them for T> because the 
probabilities given by C . 



fi = fzi ° /21 hi h 

h = fz2 ° /21 hi h 

h = hi /21 fu h 

h = f%2 ° /21 fl2 h 



fll ° /21 /31 

fll ° /21 /32 

/l2 /21 /31 

fl2 ° /21 /32 



The update functions are the following: 

/i(x) = (l,x 2 , x 2 ) / 2 (x) 
/3(x) = (1, x 2 , x 2 + x 3 ) / 4 (x) 
/ 5 (x) = (1 + xi, (x! + l)x 2 , (a?i + l)x 2 ) / 6 (x) 

/ 7 (x) = (1 + Xi, (Xi + 1)X 2 , (x x + l)x 2 + X 3 ) / 8 (x) 

(4) The probabilities that we assign are: pi = C(/i) = .18;p 2 = C(/ 2 ) = .12;p 3 = 
C(f 3 ) = .18; p 4 = C(/ 4 ) = .12;j?5 = C(/ 5 ) = .12; pe = C(/ 6 ) = .08;p 7 = c(/ 7 ) = 
.12; j9g — c (/s) = -08. In order to study the state space, it is convenient to determine 
the transition matrix of the system. The transition matrix of T> is T, 



(l,X!X 2 ,XiX 2 ) 

(1,XiX 2 ,XiX 2 +x 3 ) 

(1 + X 1? XiX 2 ,XiX 2 ) 

(1 + Xi,XiX 2 ,XiX 2 +x 3 ) 



T 
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2.1.1. Example. We notice that there are several PSS that we can construct with 
the same initial data of functions and permutations, but with different set of prob- 
abilities, that is, sub-PSS of V see 13.11 For example if S = {fi, / 2 , / 3 , f^} and 
C = M/i) = .355,d(/ 2 ) = .2IMC/3) = .12,d(/ 4 ) = -314)}, then 

B=(T;F 1 ,F 2 ,F 3 ;a 1 ,a 2 ;S; C). 
3. HOMOMORPHISMS OF PROBABILISTIC SEQUENTIAL SYSTEMS 

One objective in this section is to show that: homomorphism of PSS is a natural 
extension of homomorphism of SDS . We remark, in the introductory section [H the 
definition of SDS . In the definition of a homomorphism of PSS we establish condi- 
tions to connect the support graphs, the State Space and the assigned probabilities. 



3.1. Definition of homomorphisms of PSS. Let T>\ = (V, (-Fi)| =1 (a? 
and V 2 = (A, (Gi)£i =m , (f3 k )k, S 2 , C 2 ) be two PSS over a finite field K. 
A homomorphism from T>\ to T> 2 is a pair of functions (<f), h^) where: 

(13.11 1) <p : A — > T is a graph morphism, and 

(0 6 : fC^ip) -> K, b \/b G A), 



ft, co 
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.12=p5 



pl = 


1 8; p2 = 


12; p3 = 


% 8 T P 4 = 


.12, 


p5 = 


A2;p6 = 


.08; pi = 


.12; pS s= 


.08 




.434=d3+d4 



<fl - .355;^2 = .21 1;^3 = .12^4 - .314] 



Figure 1. Above the State Space of V, and below the State Space of B 



is a family of maps in the category of Set. They induce the adjoint function 
h , see (OD (I). 

(13.11 2) The induced adjoint map h<f> : JC n — > JC m is a map such that for all update 
function / in Si there exists an update function g E S2 such that (0, h$) is an 
SDS-homomorphism from (T, (/ : /C n -> /C n ), a*) to (A, : /C m -> JC m ),p k ). 
That is, the diagrams (Dl), and (D2) commute for all / and the selected g. 

(13.11 3) (The e-condition) For a fixed real number < e < 1, the map satisfies 
the following: 

max U;V \c fi (u,v) - dg^h^u),^^))] < e 
for all / in Si, and its selected g in 5*2 by condition (13.11 2). and u, v G /C n , 
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We will say that an homomorphism (0, h^) from T>i to T> 2 is an isomorphism if </> 
and h$ are bijective functions, and d{h^{u), h</>(g(u)) = c(u,f(u)) for all u, in /C n , 
and all /. We denote it by V l = V 2 . 

3.2. Theorem. In the definition of homomorphism the condition ( Iff. 11 3) is a con- 
sequence of the condition ( Iff. II 1), and $3.11 2). 

Proof. Suppose (0, h^) satisfies (13.11 1). (13.11 2). and max U}V \cf.(u, v)— d g .{hj,(u), h${vj)\ > 
1, then for some (u, v) we have one of the following cases 

(Case 1) Cf^Ujv) = 1 and dg^h^ujjh^v)) = 0, it is impossible by (13.11 2). 

(Case 2) dg^h^u), h^{v)) = 1, and c^iu.v) = 0, it is impossible because there ex- 
ists at least other element w G fC n such that Cf^Ujiu) ^ 0, and so on 
d g .(h^(u),h^(w)) 7^ 0. Because the sum of probabilities of all arrow going 
up of h<f,(u) is equal 1, then dg.^h^u), h^(v)) < 1. 

Therefore the third condition (13.11 3) holds, and always e exists. □ 

3.3. Theorem. If (<ft, h^) : T>\ — > T> 2 is a PSS-homomorphism, and Tj denote the 
transition matrix ofT>i, and the entry (u,v) ofT\ is p(u,v) then: 

\{ T l)u,v ~ { T 2)h^u),h^v) \ = 0, 

for all possible u and v in JC n . 

Proof. The e-condition (13.11 3) asserts that for a fixed real number < e < 1, the 
map satisfies the following: 

max u , v \c fi (u,v) - d gj ((f)(u), <f>(v))\ < e 

for all / in S±, and its selected g in S 2 by condition ( 13.11 2). and u, v G /C". 

If we have a function /j going from u to v — fi(u) in /C™, then there exists a 
function gj going from h^iu) to h^v), so gjih^iu)) = h^{fi{u)). Then for m = 2, 
and by the Chapman-Kolmogorov equation [161 E], we have the following: 

\c f f(u,fi(u)) - dg^h^u), g^ih^u)))] = 

\cfi(% fi(u))c fi (fi(u), f?(u)) - dg^h^u), 0i(/ty( u ))Ki(&'(M u ))> g 2 (h^{u)))\ = 
\c fi (u, fi{u))c fi (fi(u), ff{u)) - dg^h^u), h^fiiu^dg^h^fiiu)), /ty(/?(u)))| < 
<\c } mu)jKu))\e+\dg j {h 4> {u)Mm))\e<2e. 
By condition (2) in definition of homomorphism. Then we proved that \cj2(u, ff{uj) — 
d g 2((j)(u), g 2 (4>{u)))\ < 2e. Using mathematical induction over m, we can conclude 
that 

01) maxKc^K/rH)-^^),^" 1 ^)))! <me 

Let 5*0 be the set of function in S 2 associated to a function in S±. If u, v G /C n , and 
we denote by p 2 (u,v) = Y.f? c ff( u ' v ) and P2(M U )>M V )) = J2 g 2 d g ?(h<f>(u), h^v)), 
then condition (13.11 3) implies that \p 2 (u, v) — p^h^u), h<f,(v))\ < Ike + 6 2 , where k 
is the maximum number of functions ff going from one state to another in the state 
space JC n . Denoting 5 2 = Ylg 2 ^^ dgtih^u), h^(v))we know that 5 2 < 1. So, if 
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denote the transition matrix of T>i, and the entry (u, v) of {Ti) 2 is p2(u,v) then the 
condition (13. 11 3) implies that: 

(0312) max u ,„ - (Ta)^,^,! < 2ke + 5 2 

for all possible u and t> in K, n . 

Using equations (13.31 1). and (13.31 2). and denoting 5 m = Ylg^sjldg^h^iu), h^v))]" 1 
we conclude that 

J33I3) max M ,„ | (TO™ - (T 2 )™ {u)Mv) \ < mke + 5™ 

for all possible u and v in K. n . 

So, for all real number < e' < 1 there exists m G N, such that, 

— (^2)^ > ( u ) ) / l ^(„)| < e\ 
for all natural number m' > m, and for all possible u, v G /C n . 

In fact, using notation of equation (13.31 2). we have e m ' <C e' m , and this implies 
(m'k)e m ' < (mk)e m . Similarly 5 m ' < 5 m , So, selecting m such that (mk)e m + 5 m < e', 
we obtain 

| (TO™; - (Ta)^),^)! < (m'k)e m ' + 5 m ' < (mk)e m + 5 m < e', 

where k is the maximum number of functions going from one state to another in the 
state space of the power w! of the functions . Therefore 

jisJTC, -(^ H , M ,)l = 0, 

for all possible u and v in K n , and the theorem holds □ 

3.4. Definition of e-equivalent. If <ft and h$ are bijective functions, and the con- 
dition (13.11 3) holds, but the probabilities are not equal, we will say that X> 1; and 
T> 2 are e-equivalent, and we write T>i ~ e V 2 . So, V 1 , and V 2 are e-equivalent if 
there exist (h<f,,<f)), and ((h7 , 0" 1 ), such that for all f E Si and (7 G S^, we have 

/ = hj 1 o g o 

3.5. Image of ft^. Consider d(g) defined as follows: if S^ = {gi, . . . ,g s } is the set 
of functions in S 2 selected in (13.11 2) for the map then the new probability of g^ is 

./ \ d (9i) 

d{9i) = ^ "77 7- 

The set of functions S$ in Q, together with the new probabilities defined above, form 
a new PSS, that we will call Image of = Im^h^). So, the graph A, the update 
functions determine de local functions associate to each vertex of A, similarly the 
permutations using by these functions, and finally the new probabilities assigned. 

3.6. Sub Probabilistic Sequential System. We will say an injective monomor- 
phism is a PSS-homomorphism such that is surjective and the set of functions 
4>i, for all % are injective functions, and so h^, [llj. Therefore, we will say that 
a PSS Qx is sub Probabilistic Sequential System of Ts if there exists an injective 
monomorphism from Qx to Ts- 
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3.7. Definitions. 

(157711) Let T = (T, (F;)™ =1 , (a j ) jeJ , C) be a PSS. The pair of functions J = 
(idr,id)c-n) is the identity homomorphism, and it is an example of an isomorphism. 

(13.71 2) An homomorphism (0 l5 h^) of PSS is an injective monomorphism if is 
surjective and /i^ is injective, for example see (HJ3). Similarly we will say that an 
homomorphism is a surjective epimorphism if is injective and is surjective, 
for a complete description of the properties of this class of monomorphism and 
epimorphism see Section 7 in [TT] . 

3.8. Simulation of PSS. We consider that the PSS Q is simulated by JF if there 
exist a injective monomorphisms (01,/i^J : F — > G or a surjective epimorphism 
(02, V) : £? -» -T 7 - 



4. Examples of Simulation and Homomorphisms of PSS. 

In this section we give several examples of PSS-homomorphism, and simulation. 
In the second example we show how the condition (13.11 2) is verified under the 
supposition that a function is defined. So, we have two examples in (jH2), one with 
the natural inclusion, and the second with the only possibility of a surjective 
map. In the last example we have a complete description of two PSS, where we have 
only one permutation and two or less functions for each vertices in the graph. In 
particular this homomorphism is an injective monomorphism, so is an example of 
simulation too. 

(Hll) For the PSS, in the examples 12.11 and 12.1.11 we now define the natural 
inclusion (Idr, i) : B —> T>. It is clear that the inclusion satisfies the two first 
condition to be a homomorphism, and the third one is a simple consequence of the 
Theorem 13. 2[ In fact 

1 
0.566 0.434 
0.525 0.475 
0.211 0.314 0.12 0.355 
1 
0.566 0.434 
1 
0.08 0.434 0.566 / 










\o 



and \(Tjg) U)V — (T) t ( u ) )t ( v )| < 0.4. Here V is simulated by B 
(IH 2) Consider the two graphs below 



2 — 3 



1 — 4 



and A 



Suppose that the functions associated to the vertices are the families {/i, f 2 , /3, /k}, 
for T and {gi, g 2 , #3} for A. The permutations are a>i = (4 3 2 1), « 2 = (4 1 3 2) 
and A = (3 2 1), ft = (1 3 2), so, S = {/ = / 4 o / 3 o / 2 o / l5 / = / 4 o f x o / 3 o f 2 }, 
and X = {g = g 3 o g 2 o g = g x o g 3 o g 2 }. Then, we have constructed two PSS, 



10 



MARIA A. AVINO-DIAZ 



each one with two permutations and only one function associated to each vertex in 
the graph; denoted by: 

V s = {T]fi,f2,f3,faa 1 ,a 2 ]S\C) and B x = (A; g 1 , g 2 , g 3 ; (3 l , f3 2 ; X; D). 

Case (a) We assume that there exists a homomorphism (0, h^) from V s to Bx, 
with the graph morphism : A — > F is given by (f)(1) = 1, (f)(2) = 2, 0(3) = 3. 
Suppose the functions 

(0 6 : /C^( 6 ) -> /C 6 , Vb G A), 
are giving, and the adjoint function 

hf : Z/ -> Z p 3 , ^(^1,^2,^3,^4) = (0i(^i), 02(^2), 03(^3)) 

is defined too. If (h^,4>) is an homomorphism, which satisfies the definition (13.11) . 
then the following diagrams commute: 

Z4 h v y 4 /3 , y 4 /2 v y 4 h t y 4 7 4 / rj 4 

p ' Zip ► Zip ► Zip > Zjp Zjp ► Zip 

^ J. ^ ]• ^ J. h^l , ^ [ hj, I 

V 3 H 3 93 7 3 92 7 3 91 y 3 y 3 9 , y 3 

Zjp > Zi p > Zjp ► Zi p ► Zjp Zip y Zip 

Z4 h , y 4 h s ry 4 /3 ^ y 4 /2 v y 4 7 4 I y 4 

p ► Zip ► Zip y Zjp y Zjp Zi p > Zip 

^ J, ^ J. ^ J, ^ i h<j, i , h<j> I h<t, i ■ 

y 3 Jd v 3 _£i y 3 _£3 y 3 _S2 7 3 3 9 3 

Zjp — y Zip — y Zip — y Zi p — y Zi p z p — > Zj p 

Case (b) Consider now the map (f> : Y — ► A, defined by 0(1) = 1, 0(2) = 2, 
0(3) = 3, and 0(4) = 1. If there exists an homomorphism (0, h$) : -Bx — > ^5 that 
satisfies (13.11 2. then 

rj 3 93 7 3 92 7 3 91 7 3 rj 3 9 v3 

Zjp — > Zjp — > Zjp — y Zj 2 Zjp — y Zi p 

hj, i i h^l ^ i 

/ 4 /4°/3 7 4 h rj 4 /l y 3 7 4 / y 4 

Zjp — y Zj v — y Zj v — y Zi p Zj p — > Zi p 

y 3 91 y 3 93 y 3 92 y 3 y 3 g y 3 

Zjp — y Zip — y Zip — y Z12 Zj p — > Zj p 

h^i h<j, i h<f> i ^|,^| h^i 

/ 4 U°fl rj 4 h rj 4 h 7 3 „ 4 / ry A 

zjp ► zjp > Lip ► zjp Zp y zip 

flH3) We now construct a PSS-homomorphism from F s = (I\ (Fi) 3 , j3, S, C) to 
Qx = (A, (Gi)4, a, X, D), with the property that is surjective and the functions 
0i are injective, that we call a injective monomorphism. The PSS Ts has a support 
graph T with three vertices, and the PSS Qx has a support graph A with four 
vertices 

3 2 — 4 

r I a \ 1 , 

1 — 2 1 — 3 

The homomorphism (h^,, 0) : JF 5 — > has the contravariant graph morphism 
: A — y T, defined by the arrows of graphs, as follows 0(1) = 1, 0(2) = 0(3) = 2, 
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and 0(4) = 3. The family of functions 0« : — > JC(i), 0i(xi) = X\\ 4> 2 {x 2 ) = 
X2'-, 03(^2) = x 2 ; 04(^4) = x 4 . The adjoint function is 

: Z 2 3 -> Z 2 4 , with h ( / ) (x 1 ,X2,x 3 ) = (0i(xi), 2 (x 2 ), 3 (x 2 ), 4 (x 4 )) = (xi,x 2 , x 2 , x 3 ). 
The first condition in the definition 13.11 holds. 

The PSS Ts = (T; (Fi) 3 ; /3; S"; C), with data of functions Fi = {/ n = Id, /i2(^i, ^2, ^3) = 
(l,x 2 ,x 3 )},F 2 = {f 2 i{xi,x 2 ,x 3 ) = (xx,x 2 ,x 3 )}, andF 3 = {/ 3 i(xi, ^2, x 3 ) = (^1,^2,^2^3)}, 
one permutation or schedule /3 = ( 1 2 3 ); and probabilities C = {c(/) = .5168, c(g) = 
.4832}, so S = {f,g}; where the two update functions are 

/ = hi /21 /31, f(xi, x 2 , x 3 ) = (xi, x 2 , x 2 x^); 

and g = f 12 o f 21 o / 31 , g(x h x 2 , x 3 ) = (1, x 2 , x 2 ~xt). 

The PSS Qx = (A; (Gi) 4 ; a; X; £>) is a PSS, with the following data: the families 
of functions: G x = {gn,gi 2 }; G 2 = {g 2 i,g 22 }, G 3 = {g 3 i,g 32 }; and G 4 = {#4}- 

gu{xi,x 2 , x 3 , x 4 ) = (1, x 2 , x 3 , x 4 ) 

^21(^1,^2,^3,^4) = (xi, l,x 3 ,x 4 ) g 12 = Id = g 22 

g3i{xi,x 2 ,x 3 ,x 4 ) = (xi,x 2 ,x 1 x 2 ,x 4 ) g a a(x 1 ,X2 1 x 3 ,x i ) = (x lt x 2 ,x 2 , x 4 ) ' 

<7 4 l(xi, X 2 , X 3 , X 4 ) = (xi, X 2 , X 3 , X 2 X 4 ) 

One schedule a = (l 2 3 4), the eight possible update functions, and its 
probabilities D = {d{g) = 0, d(f) = 0, d(f) = 0, d(g) = 0, d(f) = .00252, d(g) = 
.08321, d(f) = .51428, d(^) = .39999} whose determine X = {f,g,f,g} probabili- 
ties, are the following: 

1) g(xi,x 2 ,x 3 ,x 4 ) = (gu o g 2 i o g 31 o g 4 )(x 1 ,x 2 ,x 3 ,x 4 ) = (1, 1, xix 2 , x 2 xl), 

2) /(xi,x 2 ,x 3 ,x 4 ) = {gi 2 og 21 o g 32 o g 4 )(xx,x 2 ,x 3 ,x 4 ) = (xi, 1, x 2 , x 2 xt), 

3) /(xi,x 2 ,x 3 ,x 4 ) = (g 12 o g 21 o g 31 o # 4 )(xi, x 2 , x 3 , x 4 ) = (xi, 1, xix 2 , x 2 x7), 

4) ^(xi,x 2 ,x 3 ,x 4 ) = (g n og 21 o g 32 o g 4 )(xi, x 2 , x 3 , x 4 ) = (1, 1, x 2 , x 2 xl), 

5) /(xi,x 2 ,x 3 ,x 4 ) = (5-12 o g 22 og 31 o g 4 ) (xi, x 2 , x 3 , x 4 ) = (x u x 2 , x x x 2 , x 2 xt), 

6) ^(Xi,X 2 ,X 3 ,X 4 ) = (5-11 O g 22 O g 31 0fi ( 4 )(xi,X2,X 3 ,X 4 ) = (l,X 2 ,XiX2,X 2 X4"), 

7) /(xi,x 2 ,x 3 ,x 4 ) = (5-12 o g 22 o g 32 og 4 ) (x i,x 2 , x 3 , x 4 ) = (xi, x 2 , x 2 , x 2 x7), 

8) ^(xi,x 2 ,x 3 ,x 4 ) = (5-11 o g 21 o g 32 o g 4 )( Xl ,x 2 ,x 3 ,x 4 ) = (1, x 2) x 2 , x 2 x7), 

We claim (0, h^) : JF 5 — > Q x is a homomorphism. We will prove that the following 
diagrams commute. 

Z3 — s ry 3 V 3 — s ry 3 

2 > lx2 ^2 > A 2 

h^i i , and [ J, hj, , 

Z 2 4 M Z 2 4 Z 2 4 -±> Z 2 4 

In fact, 

{h<f,0 f) (Xi,X 2 , X 3 ) = ^(xi,X 2 ,X 2 Xi) = (Xi,X2,X2,X 2 X3~) = (/o ^)(xa,x 2 ,x 3 ). 

On the other hand, 

(/i^og)(xi,x 2 ,x 3 ) = h(l,x 2 ,x 2 x^) = (l,x 2 ,x 2 ,x 2 xi) = o /i )(xi,x 2 ,x 3 ). 
We verify the composition of functions as follows 
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Figure 2. State Spaces of T and T' . Transition Matrices: T x , and Th 
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To describe the connection between the state spaces, see Figure 2. 



5. The category PSS 



In this section, we prove that the PSS with the homomorphisms form a category, 
that we denote by PSS. In Theorem 15.31 we prove that the category of Sequential 
Dynamical Systems SDS is a full subcategory of PSS. 

Theorem and Definition 5.1. Let 



be two PSS-homomorphisms. Then the composition Tt = TL2 °TL\ ■ T ' — > C is defined 
as follows: Tt = (0i o 2 , h<j, 2 o h^). 

The map Tt : T — > C is a PSS-homomorphism. 

Proof. The composite = 02 ° <f>\ '■ A —>■ T of two graph morphisms is obviously 
again a graph morphism. The composite = o is again a digraph morphism 
which satisfies the conditions (13.11 1). and (13.11 2). In fact, using the Proposition and 
Definition 2.7 in [TT], we can check these condition and conclude that the third 
condition holds, too by Theorem 12.1.11 So, (0, h^) is again a PSS-homomorphism 



Theorem 5.2. The Probability Sequential Systems together with the homomor- 
phisms of PSS form the category PSS. 



Theorem 5.3. The SDS together with the morphisms defined in VJ.1\ form a full 
subcategory of the category PSS. 



In paper [TT], the authors proved that the category SDS has finite product, simi- 
larly the category PSS has finite product too, so the following theorem for two PSS 
holds [12] 

Theorem 5.4. Let T>s, and Qx be two PSS over the finite field K,. For all e- 

homomorphisms : C — > T>s and (4>2,h^ 2 ) : C — > Qx, then there exists a 

morphism (0, h^) : £ T>s x Qx such that the following diagram commutes 



Hi = (0i, hfa) : T ^ Q and H 2 = (02, hfa) : Q — ► £ 



□ 



Proof. Easily follows from Definition and Theorem 15.11 



□ 



Proof. It is trivial. 



□ 




(<f>l, h <hi ) ^ (<j>2,hd>o) 
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5.1. The functor T : PSS -> Ab. Let V = (T, (F;)^, ( aj ) j} S,C) be a PSS 
over the finite field K. Let F = {/ : /C n -> /C n } = SOS be the set of all update 
functions that we can construct with the local functions and the permutation in 
J-'s- Let us consider, the free abelian group generated by F, that we denote by 
(F), then (F) = (S) © (S). We can notice that we are working over a finite field 
with characteristic the prime number p, so pf = for all / £ S. We will take the 
quotients of these groups by p, that is (F) j (pF) = (S) / (pS) © (S) / (pS), and these 
groups are finite, [H El E] ■ Denoting (A) / (pA) = (A) p , for an abelian group A, we 
rewrite the above relation by (F) p = (S) P Q)(S) P and, there exists a covariant functor 
from the category PSS to the category of small abelian groups with morphism of 
such, Ab, defined as follows. 

T : PSS -> Ab 

1. The object function is defined by T(Ts) = {S) p . 

2. The arrow function which assigns to each homomorphism (0, h</,) : {J-'s) (Gx) i n 
the category PSS an homomorphism of abelian groups T(</>, h^) = : (S) p — ► 
which is defined in a natural way, because h$ o (^j g5 af) = (X^es a 5') ° ^ w bere 
a G Z p , and C X, then 

T is a functor , in fact, 

1 T{l rs ) = l {s)p . 

2 T((0 2 , /l^) O /l^)) = T((0 2 , hfo)) O T((0 1; /l^)) 

The functor T gives the possibility to work with PSS using group theory, for example, 
because (F) p = (S) p © (S) p , and 



We assign probabilities to the set S in some way, and we consider that all possible 
different assignations are e-isomorphic. 

Definition 5.5 ((Complement of a PSS)). A complement of the PSST>s is the PSS 
T>-§, and all of the complement are e-equivalent, so we can select a distribution of 
probabilities for the complement having in account particular applications. 

We use the definition of complement in order to define a decomposition of a PSS 
in two sub PSS, only looking the set of functions. One of the mean problem in 
modeling dynamical systems is the computational aspect of the number of functions 
and the calculation of steady states in the state space, in particular the reduction 
of number of functions is one of the most important problem to solve for determine 
which part of the network state space could be simplify. 
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